NEW RESULTS ON DIRECT CP VIOLATION IN CHARGED KAON DECAYS 

BY NA48/2 



G. LAM ANN A° 
Department of Physics & INFN Pisa, 3 Largo Pontecorvo, 
56127 Pisa, Italy 



The NA48/2 result, based on the data collected during the 2003 run, on direct CP vi- 
olation in 7r^7r''7r'' decay is presented. The main goal of the experiment is to 
reach a sensitivity at level of 10~* in the measurement of the charge asymmetry parame- 
ter Ag — {g'^ — g^)/{g^ + where g is the linear slope of the Dalitz plot in the K ^ Stv 
decay. Thanks to the simultaneous collection of the two kaon charges and to the high reso- 
lution of the main sub-detectors, the systematics uncertainties are kept under the statistical 
error level. The experimental procedure, the analysis technique and the main systematics are 
discussed to present the final result 



(1.8 ± 2.6) X 10" 



This result based on more than 45 • 10^ events, correspondig to one half of the whole two 
year data taking, is about an order of magnitude more precise with respect to the previous 
measurement. 



1 Introduction 

In spite of the fact that more than 40 years have passed since the discovery of the CP violation^, 
the fuh understanding of this phenomenon is stih far to be reached. After the discovery of direct 
CP violation in the neutral kaon decay into two pions in the first part of the 90s, bv the NA31 
experiment 1^, and the clear confirmation after few years, by NA48l^ and KTEvEl, the only 
other environment in which this effect was detected is the neutral B systemEI. Qn the other 
hand, the study of the tiny effects due to the violation of this symmetry in all the systems where 
it's possible, represents an important window on the contribution of new physics beyond the 
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Standard Model (SM). In fact new effects could appear in a relevant way in the heavy quark 
loops which are at the core of the mechanism allowing CP violation in the decay. In the kaon 
sector the most promising places, besides s' je^ where this kind of contributions could play some 
role are the rates of GIM suppressed rare decays^ and the charge asymmetry between charged 
kaons. In particular in IT ^ Svr this asymmetry could give a strong qualitative indication 
of the validity of the CKM description of the direct CP violation or the existence of possible 
sources outside of this paradigm. Quantitative tests or predictions are, at the moment, very 
difficult from a theoretical point of view due to the complexity of the calculation involving 
non-perturbative hadronic effects. As soon as the intense theoretical efforts to improve the 
understanding of this sector of the particle interactions will produce tools to make reliable 
predictions, the experimental measurements could also be used as strong quantitative test. 

From an experimental point of view, the easiest way to look for any possible differences 
between K'^ and decay into three pions is to compare the shape of the Dalitz plot distribution 
instead of the decay rates. Exploiting the small Q value of the — > Svr decay it's possible to 
use a polynomial expansion of the matrix element*: 

\M{u,v)\'^ 1 + gu + hv? + kv'^ + ... , (1) 

where g,h,k are called linear and quadratic slope parameters and where the Lorentz invariant 
Dalitz plot parameters!^ u and v are related to the energy {E*) of "odd" and "even" pion (i.e. 
the unpaired pion and the same sign two pions, respectively) in the center of mass, through the 
definition: 

(•S3 -So) „ ruK/'i-El (s2-si) „ El-E^ 
u = ^ = 2mA- n , V = K — = 2mK— — ■ (2) 

where Sj = {px — Pi)'^, being px and pi the kaon and the pions (with i = 1, 2, 3 and i = 3 for 
the "odd" pion) four-momenta, sq = (sj + S2 + S3)/3, and uik, rriT^ the kaon and charged pion 
mass, respectively. The present valuel^ for the linear slope ^ appearing in are, for the two 
possible decay mode: 

^(vr^^vrV) = 0.638 ± 0.020 ^(vr^^Tr+vr") = -0.2154 ± 0.0035 (3) 

and (|/i|, |A;| <^ \g\). A difference between and g~ , the linear slopes in the and K~ decay, 
represents a signal of direct CP violation. This can be formalized defining the CP violating 
parameter Ag as: 

""'-g^ + g- 2g ' 

where Ag = g^—g~ and g is the average between g^ and g^ . In the SM framework the prediction 
of the Ag value is very difficult due to hadronic effects contribution, but all the calculations 1^ 
are in the range between 10~^ and few 10~^, both in the neutral (the — > Tr^Tr^Tr*^) and in 
the charged (the — > 7r^7r"'"7r~) decay mode. Calculations involving models beyond the SM 
allow for an enhancement up to the level of few 10~^. Several experiments in the past have 
searched for the asymmetry in "neutral" and "charged" mode as summarized in tab^ • The 
sensitivity reached so far is at level of few 10~^, dominated by the systematic uncertainties. The 
main goal of the NA48/2 experiment is to reach the sensitivity of 10~^ both in "neutral" and in 
"charged" mode, covering the gap existing between the experimental results and the "beyond 
the SM" theory predictions. 



Due to the strong interaction process ■k'^tt tv tt this simple parametrization isn t fully adequate to 
describe the matrix element"-! Nevertheless this expansion is adopted here for the CP violation measurement 
The linear slope for the v variable is forbidden for symmetry reasons. 



Asymmetry 


# 0/ events 


Experiment 


A^g = (19 ± 125) • 10-4 
^0 = (2 ± 19) • 10-4 


115K 
620K 


CERN PS(1975)III1 
Protvino IHEP (2005)1^ 


A'g = (-70 ± 53) • 10-4 
= (22 ± 15 ± 37) • 10-4 


3.2M 

54M 


BNL AGS (1970/^ 

HyperCP (2000) prelim.^'^ 



Table 1: Summary of the experimental situation both in "neutral" (Ag) and in "charged" (Ag) mode, before the 

NA48/2 results 



The NA48/2 experiment (basically an upgrade^^ of the NA48 experiment, devoted to the 
direct CP violation measurement in the system) will be briefly described, the analysis of the 
charged asymmetry in vr^vr^Tr*^ will be discussed and the final result, based on about a 

half of the whole NA48/2 data sample, will be presented. 



2 NA48/2: beam line and experimental apparatus 

A novel high intensity beam line was designed in order to allow for the simultaneous collections 
of and K~ decays. The possibility to collect at the same time, with the same detector, the 
X — > 37r decay of both charges, represents a fundamental point in the reduction of the systematic 
error. The beam transport line is schematically sketched in figlHl The hadron beam, of both 
charges, is produced by the SPS 400 GeV/c protons (~ 7-10^^ protons per pulse) impinging on a 
40 cm long and 2 mm diameter beryllium target with an angle of zero degrees. A magnetic device, 
called first achromat, selects the momentum of the beam in the range (60 it 3) GeV/c, splitting 
the two charges. After being recombined the beams are focused by a quadruplet of quadrupoles. 
The second achromat system houses a spectrometer for the beam, called K ABES to measure 
the particles' momentum with a resolution of 1%, useful in particular for rare decays studies. The 
beam, recombined again along the beam axis, contains ~ 6.4 • 10^ particles per 4.8 s spill (about 
12 times more pions than kaons) with kaon charge ratio K'^/K^ 1.8, which is irrelevant for 
the measurement . In the 114 m long decay region the two beams, both ~5 mm wide (RMS) are 
superimp ose d at level of 1 mm. The central detector is based on the old NA48 detector described 
elsewhere^l. The spectrometer magnet was operated in order to give P±kick = 120MeV/c and 
the resolution in momentum (GeV/c) is o'p/p = 1.0% © 0.044% • p. The LKr electromagnetic 
calorimeter is employed to detect the photons from vr'^ decay with an energy resolution (GeV) 
of aE/E = 3.2%/VE © 9%/£' © 0.42%. Thanks to these performances the resolution on the 
reconstructed kaon mass is good enough {0.9MeV/c'^ in the "neutral" mode, l.JMeV/c^ in the 
"charged" mode) for a precise calibration of the detector instabilities. 

A two level trigger is employed to reduce the rate of data collected. A hardware level trigger 
(LI) using informations from the fast hodoscope counter and from a dedicated LKr read out 
system is followed by a L2 section based on processors for fast reconstruction of kinematical 
quantities. In particular for the — > vr^vr'^Tr'^ online selection, the missing mass of the charged 
pion is required to be far from the vr'' mass, in order to reduce the contribution of the more 
frequent tt^tt" decay. The final trigger rate is 10 KHz. 



3 Strategy of the measurement 

The method to extract Ag {Ag for the "neutral" mode) is based on the comparison between the 
u projections of the Dalitz plot distribution. The ratio between the two distributions is well 
described by: 

p. ^ Nk+ 1 + (ff + Ag)u + hu^ 

Riu) = — oc ; K . (5) 
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Figure 1: The NA48/2 beam line and detector. Not to scale. 



The presence of magnetic fields botli in the beam sector (achromat) and in the detector (spec- 
trometer magnet) introduces an intrinsic charge dependent acceptance of the apparatus. In 
order to equahze this asymmetry the main magnetic field are frequently reversed during the 
data taking (the spectrometer magnet on a daily basis, while the achromat system on a weekly 
basis). For each achromat polarity, it's possible to define two ratios in which the same side (the 
two sides are called Jura (J) and Saleve (S)) of the spectrometer is involved: 

Rj{n) = -^, Rs{u) = -^ . (6) 

In each of these ratios the acceptance asymmetry due to the spectrometer is cancelled. The full 
cancellation of the detector asymmetry is obtained exploiting the quadruple ratio: 

R{u) = RusRujRnsRnj -R{ '^ |^ + ^f, + , (7) 

1 + gu + hw^ 

where U and D stands for the path followed by the in the achromat (up and down) and 
R is an inessential normalization constant. This method is independent from the relative size 
of the four samples collected with different fields configuration and from the and K~ flux 
difference. The method allows for the cancellation of: 

• local detector bias (left-right asymmetry), thanks to the fact that each single ratio in Q, 
is defined in the same side of the detector; 

• beam local biases, because in each single ratio the path of the particle through the achromat 
is the same; 

• global time variation biases, because the decays from both charges are collected at the 
same time. 

The result remains sensitive only to the time variation of the detector with a characteristic time 

smaller than the inversion period of the magnetic field. In fact each of the ratios employed in the 

((Tj) is constructed using data of subsequent days. Other systematic biases induced by effects not 

cancelled by the alternation of the magnetic field (for instance the Earth's magnetic field and any 

misalignment of the spectrometer) have to be carefully corrected. The intrinsic cancellations of 

the acceptance asymmetries in the quadruple-ratio allow to avoid the use of simulation input for 

1 8 

the measurement. Nevertheless a GEANT3 based MonteCarlo, including the time variation 
of beam and detectors, is used for systematic studies and as a cross-check of the result. 



4 Data Analysis in —>■ vr^'^Tr^vr'^ mode 

The reconstruction of the event in the "neutral" mode exploits the LKr and the spectrometer. 
In the preselection phase the tracks and the photons are selected according to quality criteria. 
The fiducial region of the detectors is chosen to avoid edge effects or miscalibrated regions. In 
particular the spectrometer's inner radial cut is chosen according to the actual beam position, 
whose displacement during the run is measured by three-track events. The decay vertex is 
recostructed from the 7 impact positions on the LKr, for each vr*^, by using the relation: 



Among all the possible 7 pairings, the two pairs with minimum vertex difference are selected 
and the final decay vertex is obtained as the arithmetic average of the two vertices. The kaon 
invariant mass is used to select the charged pion and to identify a good vr^vr'^Tr'', by 

requiring \Mi^± — M^± I < 6MeV/c^. The measurement of the charged pion momentum, 

^PDG 

slightly biased by variable DCH misalignment, is corrected by exploiting the condition Mx+ = 
Mj^- in the reconstruction of the — > 7r^7r"'"7r~ mode. The same decay mode is used to follow 
the variation of the absolute value of the spectrometer magnetic field after each periodical 
inversion, monitored online only at level of ~ 10~^. The position of the photons impact point is 
corrected to take into account the projective geometry of the calorimeter. These effects gives, 
after the correction, a negligible contribution to the systematics error. 

A potential source of systematic bias is the trigger. This was studied directly from the 
data with downscaled sample of control triggers, uncorrelated with the main trigger. In the 
"neutral" mode the LI trigger level is essentially composed by a coincidence of a signal in the 
scintillating hodoscope (Ql) and a pattern compatible with 4 photons on the LKr (NTPEAK). 
The Ql efficiency is obtained by using all the events with one track in the charged hodoscope, 
re- weighted according to the distribution of the Svr events. The inefficiency (~ 0.25%) is constant 
during the run and the upper limit of the systematic error on /S.g given by this part of the trigger 
is 0.1 • 10-^. 

The largest systematic uncertainty (1.3 • 10~^) comes from the neutral part of the LI. The 
NTPEAK inefficiency is different at the beginning (~ 0.7%) and at the end (~ 3%) of the data 
taking. The estimation of this systematics is limited by the statistics of the control sample. 
The L2 inefficiency (~ 5.7% on average during the run) is mostly (~ 70%) due to the DCH 
wires inefficiencies and was studied by using MonteCarlo simulation, including the detailed time 
variation of the spectrometer's behaviour. The systematic error assigned to Ag' is 0.4 • 10~^, 
including also other effects like data buffer overflows, algorithm inefficiency and synchronization 
effects in the trigger chain. All the sources of systematics are summarized in table 0] . 

During the 2003 data taking (one half of the whole statistics) 40.3- lO'^ K+ and 16.9- 10*5 K' 
decays in the "neutral" Svr mode have been selected for the asymmetry measurement presented 
in this paper. The final result in terms of Ag (related to A(7 through Q ) , obtained as average 
of three independent analyses that agree within few 10~^, is: 



= (1.8 ± 2.2stat. ± l■^stat.{tr^g) ± 0.8,^,,*. ± 0.2e,t.) • lO'^ = (1.8 ± 2.6) • IQ-^ , (9) 



where the external error depends on the uncertainties on the value of §{71 tt^tt^) (~ 3.1%)!^. 
This result, fully compatible with the SM predictions, is almost one order of magnitude better 
than the previous measurement by other experiments'^,'^. The ongoing analysis of the 2004 
sample will improve the statistical error of a factor of ~ 1.4, with a similar reduction in the 
systematic error. 




(8) 



Systematic effect 


Effect on Ag ■ 10"^ 


U calculation Sz fitting 


±0.4 


LKr non linearity 


±0.1 


Shower overlapping 


±0.5 


Pion decay 


±0.5 


Spectrometer Alignment & Momentum scale 


< ±0.1 


Beam Geometry 


±0.3 


Accidentals 


±0.2 


LI Trigger: Ql 


±0.1 


LI Trigger: NTPEAK 


±1.3(stat.) 


L2 Trigger 


±0.4 


Total systematic uncertainty 


±1.0 


Trigger statistical uncertainty 


±1.3 



5 Result in charged mode 

An analogous analysis in th e "c harged" mode using 1.6 • 10^ — > Tr^vr^vr^ decays, from the 
2003 sample, gives the result!^ 

A'g = (1.7 ± 2.1stat ± 1.4t„g ± lAsyst) ■ 10-^ = (1.7 ± 2.9) • 10'^ . (10) 

The statistical precision is similar with respect to the result given in ® because, even though 
the "neutral" mode is statistically disfavored due to lower branching ratio and acceptance, the 
population density of the Dalitz plot is more favourable and |5(7r^7r'^7r'^)| ~ 3 • |5f(7r^7r'''7r~)|. A 
preli mina ry study on the combined 2003 and 2004 data (3.11 • 10^^ events) leads to a compatible 
resultlSni 

^^ = (-1.3±1.5,tat±0.9t„3±1.4,j,,t)-10-^ = (-1.3±2.3)-10-4 . (11) 
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